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ABSTRACT 


The  problem  of  allocating  different  types  of  jobs  among  several 
classes  of  military  manpower  is  becoming  increasingly  more  important  as 
military  jobs  grow  more  specialized.  A  model  was  proposed  by  Richard  C. 
Grinold  which  constructs  a  personnel  inventory  by  rank  for  each  of 
several  classes  of  manpower  and  then  allocates  that  inventory  to  meet 
billet  requirements.  The  model  is  designed  for  long-range  planning 
purposes  and  produces  possible  inventory  based  on  an  optimization  scheme 
that  sets  permitted  errors  in  the  allocation.  The  thesis  presents  a 
review  and  demonstration  of  the  model  based  on  the  U.S.  Navy  officer 
corps,  a  discussion  of  implementation  considerations,  and  further  work  on 
optimization  schemes. 
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I.  INTRODUCTION 


The  problem  of  allocating  different  types  of  jobs  among  several 
classes  of  military  manpower  is  becoming  increasingly  more  important  as 
military  jobs  become  more  specialized.  A  model  was  proposed  by  Grinold 
[Ref.  1]  which  constructs  a  personnel  inventory  by  rank  for  each  of 
several  classes  of  manpower  and  then  allocates  that  inventory  to  meet 
billet  requirements.  The  model  is  designed  for  long-range  planning 
purposes  and  produces  possible  inventories  based  on  an  optimization 
scheme  that  sets  permitted  errors  in  the  allocations. 

The  next  section  reviews  the  model  and  offers  several  points  to 
consider  concerning  its  implementation  based  on  practical  problems  during 
the  preparation  of  this  thesis.  The  review  highlights  the  mathematical 
foundations  of  the  model  and  defines  the  essential  variables  and  4 

parameters.  A  complete  discussion  of  the  mathematical  concepts  is  found 
in  Grinold  [Ref.  1]. 

Section  III  presents  a  means  for  measuring  the  error  between  reauire- 
ments  and  inventories  and  then  expresses  this  measure  as  the  objective  of 
3  types  of  mathematical  programs:  linear  (goal)  programming,  elastic 
programming,  and  quadratic  programming.  Sections  IV,  V,  and  VI  present 
a  variety  of  possible  formulations  for  the  three  programming  models  ^om 
which  analysis  may  be  accomplished.  The  formulations  presented  incluae 
those  proposed  oy  Grinold  [Ref.  1]  plus  several  new  ones.  For  the 
quadratic  programming  models,  analytical  solutions  are  derived  -'rom  wnicn 
computer  programs  in  the  API  language  were  written. 
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Section  VII  is  a  demonstration  of  the  model  using  both  APL  programs 
and  a  linear  programming  optimizer  available  at  the  Naval  Postgraduate 
School.  The  demonstration  is  based  on  typical  data  from  the  U.S.  Navy 
officer  corps.  The  last  section  summarizes  the  results  of  the  investigation 
and  proposes  areas  for  further  studies. 
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II.  REVIEW  OF  THE  MODEL 


A.  PURPOSE 

The  purpose  of  the  model  is  to  allocate  different  types  of  jobs 
among  several  classes  of  manpower  in  a  steady  state  system.  The  term 
"allocation"  refers  to  a  method  for  distributing  jobs  among  several 
classes  of  people  or  people  among  different  types  of  jobs.  The 
allocations  are  determined  by  arrays  of  sharing  fractions  which  are 
defined  below.  These  sharing  fractions  represent  an  actual  or  possible 
management  policy. 

The  term  "steady  state"  refers  to  the  equilibrium  condition  prevailing 
in  the  model.  This  condition  arises  from  the  assumption  of  a  fixed 
longitudinal  manpower  flow  in  each  manpower  class  used  for  the  analysis. 
This  concept  is  defined  for  the  model  in  Grinold  [Ref.  1],  and  is  discussed 
in  detail  in  Grinold  and  Marshall  [Ref.  2],  and  in  Bartholomew  and  Forbes 
[Ref.  3].  Further  discussion  of  the  steady  state  conditions  is  found  in 
the  section  which  defines  the  steady  state  variables. 

The  final  item  of  interest  in  stating  the  purpose  deals  with  the 
planning  aspects  of  the  model.  Since  this  analysis  is  done  for  steady 
state  conditions,  it  can  be  used  most  effectively  to  test  current  or 
alternate  management  policies  for  their  long-term  effect  and  provide 
answers  to  the  "what  if"  policy  questions.  As  an  aggregate  planning 
tool,  the  model  is  intended  to  test  long-term  policies  that  may  be  used 
later  for  daily  operation  of  a  manpower  system.  Since  conditions  of 
steady  state  may  be  significantly  different  from  current  conditions, 


the  analyst  should  be  aware  of  the  limitations  of  the  model.  Often, 
the  relative  differences  between  competing  policies  or  the  magnitude  of 
the  effect  of  a  policy  may  be  the  point  of  interest  for  which  an  analysis 
is  being  performed.  The  model  should  be  used  with  these  considerations 
in  mind. 

The  model  is  motivated  by  a  study  of  the  U.S.  Navy  officer  corps 
and  examples  from  this  manpower  system  will  be  used  when  explaining 
concepts  as  an  aid  to  better  understanding  of  the  model.  A  description 
of  the  structure  of  the  U.S.  Navy  officer  corps  for  the  unrestricted  line 
officers  is  found  in  NAVPERS  15 197 A  [Ref.  4]. 

3.  DEFINITIONS  ANO  TERMINOLOGY 

1.  Attributes 

The  definition  and  understanding  of  variables  in  the  model  is 
complicated  by  the  large  number  of  variables  present.  The  most  important 
variables  relating  to  concepts  are  presented  in  this  section.  Other 
variables  relating  to  specific  optimization  techniques  or  used  for 
mathematical  clarity  will  be  introduced  as  the  need  for  them  arises.  A 
list  of  variables  used  in  the  model  is  found  in  Appendix  A. 

The  requirements  and  the  inventories  possess  three  common  attributes. 
They  are: 

•  MANPOWER  CLASSES--the  types  of  manpower;  indexed  by  K=l,2, . . . ,K. 

•  CAREER  STAGES--the  experience  levels  required;  indexed  by 
i  =  l,2, . . . , I . 

•  J08  CL4SSES--the  types  of  joos  to  be  filled;  indexed  by  j=l,2,...,J. 

The  clear  definition  of  attributes  is  closely  related  to  the  puroose 
of  analysis.  In  the  demonstration  wnich  follows,  the  manpower  classes 


represent  types  of  officers  (e.g.,  Pilot,  Naval  Flight  Officer);  the 
career  stages  represent  ranks  (e.g.,  Lieutenant,  Commander);  and  the  job 
classes  are  naval  officer  billet  categories  (e.g.,  1310,  1320). 

The  level  of  detail  is  also  closely  related  to  the  purpose 
of  analysis.  For  example,  the  manpower  class  of  pilot  could  be  further 
broken  down  into  pilots  by  type  of  aircraft,  but  this  imposes  increased 
data  requirements  which  may  not  be  needed  if  the  purpose  of  analysis  is 
to  study  Navy-wide  trends  rather  than  issues  specific  to  the  Navy's 
aviation  community. 

The  attributes  may  apply  to  either  requirements  or  inventories. 

In  Navy  terminology,  requirements  are  usually  stated  in  terms  of  "billets." 
Some  confusion  may  arise  since  this  is  synonymous  with  the  term  "job." 

The  reader  should  bear  in  mind  that  requirements  may  be  referred  to  as 
billets  or  jobs  which  are  distinct  from  the  attribute  of  "job  type." 
Likewise,  the  term  "inventories"  is  used  interchangeably  with  stocks, 
personnel,  and  people,  but  is  distinctly  different  from  the  term  "manpower 
types."  The  use  of  these  terms  should  be  clear  from  the  context  in  which 
they  occur. 

The  demonstration  highlights  another  problem  found  in  practical 
application.  The  types  of  people  and  types  of  jobs  are  often  expressed 
by  the  same  terminology.  For  example,  1319  is  both  an  officer  designator 
and  a  billet  designator  for  the  manpower  class  "pilot”  and  for  the  job 
class  "pilot."  To  avoid  confusion,  the  word  equivalent  "pilot"  will  be 
used  when  speaking  of  1310  officers  and  the  billet  designator  "1310"  will 
be  used  when  speaking  of  a  corresponding  type  of  job. 


Other  interpretations  of  the  model  Dased  on  current  military 
personnel  problems  are  possible.  For  example,  manpower  classes  could  be 
defined  as  enlisted  personnel  by  mental  group  category  and  the  purpose  of 
analysis  could  be  to  study  the  effects  on  billets  of  increasing  require¬ 
ments  for  the  higher  mental  groups  with  a  shrinking  manpower  pool  from 
which  to  provide  accessions.  Another  possibility  would  be  to  define 
manpower  classes  by  sex  and  job  classes  by  types  of  sea  or  shore  billets 
and  then  study  the  implications  of  increased  accession  of  women  on 
sea-shore  rotation  cycles. 

2.  Steady  State  Variables 

Let  ni  oe  the  number  of  time  periods  an  individual  must  remain 
in  the  organization  before  entering  stage  (i+1).  For  example,  naval 
officers  in  the  rank  of  lieutenant  (jg)  will  remain  in  the  Navy  four 
years  before  being  promoted  to  the  rank  of  lieutenant,  thus  n.  for 
lieutenant  (jg)  is  four.  The  total  amount  of  time  a  person  may  spend  in 
stage  i  is: 

ni  -  ni-l  <n0  =  0)- 

An  example  of  the  amount  of  time  spent  in  each  career  stage  of  the  Navy 
officer  corps  will  be  given  in  Section  VI  (Fig.  4)  for  the  demonstration. 
The  steady  state  variables  are  defined  as  follows: 

•  Survivor  fractions  by  rank  and  manpower  class--this  is  the  matrix 

^  h 

S  indexed  over  (I,K).  The  (i,k)1'  element,  sik>  represents  the 
fraction  of  type  k  people  who  remain  in  the  organization  more  than 
’  stages.  The  element  may  also  be  interpreted  as  the 
probability  that  a  type  k  person  remains  in  the  organization 
longer  than  ni  time  periods. 


t'  h 

•  Expected  stage  length  is  the  matrix  W  whose  (i,k)  element 
represents  the  amount  of  time  a  person  in  manpower  class  k  is 
expected  to  spend  in  stage  i.  W  may  also  be  defined  as  the 
product  of  the  probability  of  a  type  k  person  reaching  stage  i  and 
the  length  of  stage  i : 

wik  =  Sik<ni  -  ni-l>- 

A  more  complete  discussion  of  survivor  fractions  and  expected  stage 
lengths  may  be  found  in  Grinold  [Ref.  1]  and  Grinold  and  Marshall  [Ref.  2 

3.  Requirements 

Requirements  are  expressed  in  three  ways: 

•  Billets  by  rank  and  job  type--the  number  of  jobs  at  each  stage  to 
be  filled.  This  is  a  matrix  B  whose  elements  are  indexed  by 

•  Billets  by  rank  and  people  type--the  number  of  people  required  at 
each  stage.  This  is  not  the  normal  way  in  which  requirements  are 
stated  and  should  not  be  confused  with  personnel  inventories. 

This  is  a  matrix  P  whose  elements  are  indexed  by  (i,k). 

•  Target  al locations--the  number  of  (i,k)  people  desired  in  (i,j) 
jobs.  This  is  a  three-dimensional  array  T  with  elements  indexed 
by  (k,i,j). 

4.  inventories  and  Accessions 

Finding  optimal  inventories  and  accessions  will  normally  be  the 
qoal  of  this  model.  They  are  defined  below. 

*■  ft 

•  Accessions  are  denoted  by  a  vector  7  whose  < '  element  represents 
the  number  of  oeople  entering  manpower  class  <  eacn  year. 


Inventories  are  expressed  in  three  ways: 


•  Inventories  by  rank  and  job  type--the  matrix  X  wnose  elements  are 
indexed  by  (i,j)  defines  the  number  of  jobs  at  each  stage  actually 
filled.  This  is  not  the  normal  way  in  which  inventories  are 
stated  and  should  not  be  confused  with  billet  requirements. 

•  Inventories  by  rank  and  people  type--the  actual  number  of  available 
people  py  manpower  class  and  rank.  This  is  the  matrix  Z  whose 
elements  are  indexed  by  (i,k). 

•  Actual  al location--the  actual  number  of  (i,k)  people  who  will  fill 
( i , j )  jobs.  This  is  the  three-dimensional  array  A  with  elements 
indexed  by  (k,i,j).  The  notation  for  requirements  and  inventories 
is  summarized  in  Table  1. 


TABLE  1 

REQUIREMENTS  AND  INVENTORIES 


Job  Type 
and  Rank 

People  Type 
and  Rank 

People  Type,  Rank 
and  Job  Type 

Requirements 

bij 

Pik 

Hij 

Inventories 

Xij 

zik 

akij 

5.  Sharing  Rules 

There  are  two  rules  used  to  distribute  the  (i,j)  ioos  among  the 
( i ,«)  manpower-  types  or  vice  versa.  These  rules  are  known  as  job  sharing 
and  people  sharing. 
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•  Job  sharing  is  defined  by  a  three-dimensional  array  F  whose 

( k , i , j ) ^ h  element  gives  the  fraction  of  (i,j)  jobs  that  should 
be  performed  by  manpower  type  k. 

•  People  sharing  is  defined  by  a  three-dimensional  array  G  whose 
(k,i,j)th  element  qives  the  fraction  of  (i,k)  manpower  types  who 
should  perform  job  type  j. 


C.  DETERMINATION  OF  REQUIREMENTS 

The  target  allocation  (T),  which  is  the  number  of  type  j  jobs  to  be 
filled  by  class  k  people  of  rank  i,  is  determined  in  one  of  two  ways: 

1.  By  specifying  the  desired  billets  by  people  type  and  rank  (P)  and 
the  rule  for  sharing  these  people  among  the  different  billets  (G). 
Note  that  G  times  P  distributes  the  requirement  for  the  ktn  class 
of  people  in  rank  i  over  the  job  class  j  and  produces  the  target 
allocation  T: 

9k i j  pik  “  tkij- 

2.  By  specifying  the  billets  by  rank  and  job  type  (B)  and  the  rule 
for  sharing  these  billets  among  the  different  types  of  people  (F). 
Note  that  F  times  B  distributes  the  jcn  type  of  job  in  rank  i 
over  the  people  class  k  and  also  produces  the  target  allocation 

T: 


fki j  bij  =  bki j  * 

In  order  to  maintain  consistency,  the  following  equation  must  hold: 

9kij  pik  "  tk i j  '  fkij  bij 

The  argument  from  which  this  type  of  equation  is  derived  will  be 
referred  to  as  the  "consistency  rule."  Furthermore,  all  people  classes 
must  be  filling  some  job  and  all  job  classes  must  be  filled  by  some 
people  which  means  that  the  equation: 


(2) 


Thus,  one  may  either  start  with  F  and  3  and  calculate  P  from  equation  (2) 
and  then  6  using  equation  (1),  or  start  with  G  and  P  and  calculate  3  from 
equation  (3)  and  then  F  from  equation  (1). 

For  the  Navy  example  which  follows,  it  should  be  noted  that  the  usual 
manner  of  determining  target  allocations  will  be  to  start  with  a  given 
requirement  of  job  types  by  each  ran*  (B)  derived  from  the  force  structure. 
Since  it  is  more  common  to  think  of  people  being  assigned  to  jobs  rather 
than  jobs  to  people,  a  typical  sharing  policy  is  more  easily  understood 
in  terms  of  distributing  people  over  the  required  job  types  (G).  Thus, 
it  would  be  helpful  to  be  able  to  calculate  F  and  P  based  on  B  and  G. 
Unfortunately,  this  cannot  be  done  directly.  One  solution  to  this 
problem  is  a  least  squares  approach  using  equation  (3).  First  determine 
the  least  squares  fit  values  for  P,  and  then  solve  for  F  using  equation 
(1).  The  least  squares  objective  is: 


7  f  I  ikx 
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Notice  that  this  minimization  is  accomplished  by  minimizing  each  of 
the  following  I  expressions: 


min  £ 

P  J 


pik  '  bj)2 


for  i-1,2,.. 


I. 


For  a  particular  i,  the  above  can  be  written  in  matrix  notation: 
min  (G'P  -  3)'  (G'P  -  B). 

Thus,  the  itn  column  of  ■  *  jiven  by  the  optimal  solution: 

P  -*•  :-i  .  j  1  G  B  . 


Now  calculate  f  usin<>;  >■  vision  (1)  to  obtain: 

■1 


f*U 


gkij  ?ik  (bij} 
0 


if  bid 
if  bi . 


t  o, 

=  0. 


0.  DETERMINATION  OF  INVENTORIES 

The  actual  allocation  (A),  which  is  the  number  of  people  of  type  k  in 
rank  i  who  will  be  available  to  fill  job  type  j,  is  determined  using  the 
same  consistency  rule  invoked  for  equation  (I): 

gkij  zik  "  aki j  "  f k i j  xij 

and  the  implications  similar  to  equations  (2) 
hold.  The  equation: 


(4) 


and  (3)  above  must  also 


J  g...  *  1 
*  ykij 


implies  that 


2ik '  j  fku  x<j  • 


and  the  equation: 


imp  1 ies 


k  ^  "  1 


xij  3  l  gkij  zik  ■ 


(5) 


.6) 


1 

1 
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In  addition,  the  model  is  constrained  by  the  steady  state  requirement 
that  the  expected  inventory  of  type  k  people  in  rank  i  is  equal  to  the 
number  of  accessions  of  type  k  people  times  the  expected  duration  of  type 
k  people  in  rank  i.  In  equation  form,  this  constraint  is: 


Since  inventories  may  be  computed  from  accessions  (Y)  and  expected 
stage  lengths  (W),  the  optimization  techniques  to  determine  inventories 
found  in  the  next  sections  will  usually  have  Y  as  the  decision  variable. 
If  the  inventories  (X  or  Z)  are  the  decision  variables,  then  a  method  is 
needed  to  determine  the  corresponding  accessions.  From  equation  (7), 
note  that  for  each  y^  there  are  I  equations.  Once  again,  a  least 
squares  approach  may  be  used  to  find  the  best  fit  values  of  y^  for  this 
system  of  equations.  The  objective  is: 


min 

y 


l"ik 


This  yields  the  result: 


y,  / 1  w . , 
Jk  =  ( 7  ik 


zik)  (Vik> 


-1 


for  k=l,2, . . . ,K. 
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III.  OPTIMIZATION  TECHNIQUES 


A.  PURPOSE 

This  section  discusses  the  use  of  mathematical  programming  to  deter¬ 
mine  desired  values  for  the  variables  defined  in  Section  II.  The  degree 
to  which  these  techniques  may  be  used  depends  on  availability  of  data, 
level  of  detail  in  attribute  definition,  and  desired  accuracy  of  results. 
The  scope  of  the  problem  to  be  solved  determines  the  type  of  the  optimi¬ 
zation  technique  and  the  way  in  which  the  problem  is  formulated. 

The  purpose  of  this  section  is  to  present  a  method  by  which  the 
quality  of  an  allocation  may  be  measured.  Such  a  measure  will  be  used  as 
the  basis  for  the  objective  functions  for  three  types  of  mathematical 
programs:  linear  (goal),  elastic,  and  quadratic  programming.  Formula¬ 
tions  for  each  of  these  will  be  presented  in  subsequent  sections. 

3.  FORMULATION  OF  GOALS 

In  formulating  the  optimization,  special  care  must  be  taken  in 
determining  the  objective  to  be  achieved  and  the  decision  variable  to  be 
determined.  The  overall  objective  is  to  match  actual  allocations  to 
target  allocations  with  respect  to  a  decision  variable.  The  methods 
described  here  seek  to  do  this  by  minimizing  the  weighted  percent  differ¬ 
ence  between  inventories  and  requirements.  The  decision  variable  to  be 
determined  may  be  any  one  of  the  previously  defined  variables.  All  other 
variables  are  then  treated  as  parameters. 

L’sted  below  are  the  possible  decision  variables  and  a  brief  discussion 
of  now  they  may  affect  the  optimization. 


•  Requirements--Errors  may  be  minimized  with  either  billets  by  job 
type  and  rank  (B)  or  billets  by  people  type  ana  rank  (P)  as  the 
decision  variable.  Normally  B  is  treated  as  a  parameter  derived 
from  the  force  structure.  For  example,  if  ten  aviation  squadrons 
of  a  certain  type  are  required  by  the  force  structure,  then 
various  numbers  of  officer  pilots  of  different  ranks  will  be 
needed  to  fill  the  aviation  type  billets  found  in  those  squadrons. 
On  occasion,  it  may  be  of  interest  to  allow  the  requirements  to  be 
decision  variables  and  check  the  resulting  optimal  values  given 
fixed  inventories  and  sharing  rules. 

•  Inventories  and  Access ions--The  decision  variable  may  be  either 
the  inventories  or  accessions.  Due  to  the  steady  state  assumption, 
it  is  sufficient  to  know  the  accessions  and  the  expected  length  of 
time  in  each  career  stage  in  order  to  calculate  expected  inven¬ 
tories.  For  this  reason,  the  accessions  (Y)  will  normally  be  the 
decision  variables.  This  approach  requires  that  billets  (B  or  P), 
sharing  rules  (F  or  G),  and  expected  duration  in  each  stage  (W), 

be  treated  as  parameters.  The  objective  functions  in  Grinola 
[Ref.  I]  are  of  this  type. 

•  Expected  Duration  in  Each  Career  Stage--This  variable  may  be 

the  decision  variable  if  billets,  accessions  and  sharing  rules  are 
known.  If  certain  assumptions  are  made  about  the  conditions  from 
which  an  optimal  W  is  derived,  then  it  is  possible  to  relate  the 
optimal  W  to  a  survivor  function. 

•  Sharing  Rules--Minimizing  errors  with  F  or  3  as  the  section 
variaDle  will  result  in  the  optimal  "mixing"  policy.  If  the 
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objective  is  to  find  the  optimal  policy  based  only  on  billets, 
inventories  and  expected  stage  lengths,  then  this  approach  is 
reasonable. 

C.  MEASUREMENT  OF  ERROR 

As  stated  previously,  optimization  is  performed  to  minimize  the 
weighted  error  between  requirements  and  inventories.  The  error  for 
the  linear  models  and  quadratic  models  may  be  expressed  as  the  excess  or 
deficit  of  the  three  quantities  shown  in  Table  2. 

TABLE  2 

TYPES  OF  ERROR 

:  the  error  between  targets  and 
allocations  by  job  type,  rank, 
and  people  type. 

:  the  error  between  the  inventory 
of  (i,k)  people  and  the  corres¬ 
ponding  bi 1  lets. 

:  the  error  between  the  (i,j)  billets 
and  the  corresponding  inventory. 

The  first  error  shown  involves  all  three  attributes  and  may  be 
obtained  from  either  the  second  or  the  third  one  since  these  errors  may 
be  distributed  using  tne  job  sharing  or  people  sharing  rule: 

gkij  dik  '  ' tk i j  "  ak i j 1  '  fkij  eij 
where  dik  and  are  defined  in  "able  2. 

For  this  reason,  the  last  two  types  of  errors,  those  by  rank  and  job 
type  and  those  oy  '■ank  and  people  type,  will  be  used  in  the  following 
discussion.  In  the  linear  programing  model,  it  is  important  to  know  if 
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the  error  is  a  shortfall  or  overfill.  In  keeping  with  the  notation  in 
Lee  [Ref.  5],  shortfalls  will  be  denoted  by  the  superscript  minus  and 
overfills  by  the  superscript  plus. 

Weighting  factors  are  used  to  allow  trade-offs  between  billets  of 
differing  attributes.  The  weighting  factors  are  based  on  total  number  of 
billets  to  be  filled  and  the  relative  importance  of  Dillets.  A  shortfall 
of  twenty  billets  of  a  particular  type  will  have  a  different  impact 
depending  on  whether  there  were  100  or  1,000  billets  to  be  filled.  The 
error  values  will  have  a  common  metric  if  they  are  converted  to  percent 
errors.  Thus,  a  better  measure  of  the  error  would  be  (z-p)/p  or  (x-b)/b. 
However,  the  actual  measure  may  still  differ  depending  on  the  nature  of 
the  job  type,  rank,  or  personnel  type.  For  example,  a  five  percent 
shortfall  of  pilots  who  are  of  the  rank  captain  may  not  be  as  serious  as 
a  five  percent  shortfall  in  pilots  who  are  lieutenants.  Thus  the  percent 
shortfall  in  meeting  a  critical  target  should  be  weighted  more  heavily 
than  that  of  a  less  critical  target.  In  order  to  handle  this  problem, 
additional  parameters  indicating  the  percent  trade-off  between  ranks,  job 
types,  and  people  types  must  be  introduced  to  the  model.  One  way  in 
which  this  may  be  done  is  to  specify  parameters  which  will  serve  as 
benchmarks  to  measure  the  permitted  range  of  unit  error  in  overfilling  or 
underfilling  a  requirement. 

For  example,  let  a  five  percent  unit  error  be  the  permitted  shortfall 
of  1310  (pilot)  lieutenant  jobs.  Now  for  another  job,  say  1310  captains, 
the  relative  trade-off  for  an  underfill  may  be  found  by  answering  the 
question,  "What  percent  under  target  in  the  assignment  of  1310  captain 
billets  is  as  serious  as  a  five  percent  shortfall  in  billing  1310 
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lieutenant  billets?"  In  this  way,  essential  judgments  about  the  trade-off 
between  differing  ranks  and  job  types  may  be  made.  The  same  process 
could  also  be  done  for  ranks  and  people  types  or  the  consistency  equation 
could  be  used  to  produce  the  trade-off  parameters.  The  notation  for  these 
parameters  is  shown  in  Table  3.  These  trade-offs  may  now  be  added  to  the 
penalty  function  by  multiplying  the  percent  error  by  the  inverse  of  the 
permitted  unit  error.  Thus  the  actual  weighting  factors  for  each  error 
is  the  inverse  of  the  permitted  unit  percent  error  times  the  requirements. 
Let  u  and  v  be  the  weighting  factors  for  job  types  and  people  types, 
respect i vely.  As  before,  let  the  superscript  plus  denote  weights  associ¬ 
ated  with  overfills  and  the  superscript  minus  denote  weights  associated 
with  shortfalls.  The  notation  for  these  parameters  is  in  Table  4.  The 
total  penalty  value  is  now  defined  as  the  weighting  factor  times  the 
error.  This  penalty  value  is  a  measure  of  error  that  takes  into  account 
both  the  varying  sizes  of  the  different  requirements  and  their  relative 
importance.  The  notation  for  the  penalty  measure  deriv*-:-5  cnus  far  is  in 
Table  5. 

Obtaining  valid  weights  is  of  particular  importance  to  the  model 
since  the  final  result  of  the  optimization  could  be  very  sensitive  to 
these  parameters.  There  are  several  ways  in  which  they  can  be  estimated. 
Expert  judges  could  arbitrarily  select  either  the  parameters  to  be  used 
or  the  range  of  minimal  values  for  excess  and  deficit.  For  example,  the 
requirement  for  a  particular  job  type  and  rank  may  be  100.  If  a  permitted 
range  of  error  is  95  to  110,  then  »  =  5  percent  and  9  =  10  percent.  The 
range  of  error  is  a  useful  concept  and  the  range  bounds  a  £,  p  and  £  are 
defined  in  Taole  6.  If  the  lower  unit  errors  are  greater  than  100 
percent,  then  it  would  oe  possible  to  have  a  negative  value  for  the  lower 
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TABLE  3 


PERMITTED  UNIT  PERCENT 

ERROR 

by  job  type 
and  rank 

by  people  type 
and  rank 

shortf al 1 

*U 

5ik 

overfill 

8  .  . 

1J 

*ik 

TABLE  4 

WEIGHTING  FACTORS 


by  job  type 

by  people  type 

and  rank 

and  rank 

shortf al 1 

uij  =  ^ij  bij^ 

vik  ■  (*ik  Pik> 

overf i 1 1 

Uij  =  (9ij  bij^ 

vik  *  (‘.-k  pfk> 

TABLE  5 

WEIGHTED  PENALTY  MEASURE 

by  job  type  by  people  type 

and  rank  and  rank 


shortf al 1 


overt i 1 1 


- 

-  Xij 

_  Pik 

’  Zik 

ij 

"ij 

&ij 

;ik  Pik  5  . , 

ik 

Pik 

+ 

-faij 

,,+  :ik 

*  pik 

'ij 

9  .  . 

-J 

b<j 

Vik  3ik  » 

i  < 

Pik 
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TABLE  6 


ERROR  RANGE  OF  REQUIREMENTS 


by  job  type 

by  people  type 

and  rank 

and  rank 

lower  error 

bound 

(1  ‘  *!j>  °ij  '  i-j 

(1 

-  ‘ik1  pik  “  £ik 

upper  error 

bound 

(1  ♦  9,j)  t>(J  *  bfj 

(1 

*  *ik>  ptk  =  pik 

bound.  To  avoid  this  possibility,  it  will  be  required  by  definition  that 
the  lower  bounds  be  greater  than  or  equal  to  zero.  It  is  also  possible  to 
estimate  the  permitted  unit  error  by  job  type  if  given  people  type  and 
vice  versa  using  the  consistency  rule  resulting  in  the  equations: 

£ik  9ktJ  *  pij  fkij  for  5,1 

pik  3kU  •  5fj  fkij  for  all  k.l.j. 

D.  OBJECTIVE  FUNCTIONS 
1.  General 

Let  E . j  be  the  penalty  value  function  of  errors  in  job  class  j 
and  rank  i,  and  let  0.^  be  the  penalty  value  function  of  errors  in 
people  class  k  and  rank  i.  Then  three  possible  mathematical  programming 
objective  functions  would  be: 

(1)  min  z  r  E • • ,  i.e.,  minimizing  error  by  job  types; 

i  J  J 

(2)  min  z  z  0..  ,  i.e.,  minimizing  error  by  Deople  types; 

i  k 

(3)  min  z  z  E-.  +  z  z  0.,  ,  i.e.,  minimizinq  error  by  job  tvpes 

i  j  J  i  k  1K 

and  people  types. 

Minimizing  the  third  objective  function  should  produce  a  close  match  between 
targets  and  actual  allocations  and  was  used  in  the  mathematical  programs  in 
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Grinold  [Ref.  1J.  The  use  of  objective  (3)  requires  more  data  input 
for  parameters  than  either  objectives  (1)  or  (2).  Objective  functions 
based  on  each  of  the  three  representations  will  be  used  in  the  following 
sections . 

In  order  to  shorten  mathematical  notation  in  the  following  sections 
on  specific  objective  functions,  it  will  be  assumed  that  the  decision 
variable  of  interest  is  the  inventory  X  and  the  objective  function  is 
type  (1)  from  above.  Also,  the  subscripts  will  be  dropped  for  clarity. 

2.  Linear  (Goal)  Programming 

One  objective  function  that  measures  the  weighted  percent  errors 
is  a  piecewise  linear  convex  function  that  has  a  value  of  zero  if  the 
actual  allocation  meets  the  target,  and  increases  at  a  rate  inversely 
proportional  to  the  permitted  error  in  meeting  that  target.  This  means 
the  function  is  unity  when  the  decision  variable  assumes  either  the  lower 
or  the  upper  error  bound.  Therefore,  the  function  will  satisfy  the 
following  three  conditions: 


(i) 

E(x) 

=  0 

if 

x  =  b  , 

(H) 

E(x) 

=  1 

if 

x  =  _b  , 

(ill) 

E(x) 

=  1 

if 

x  =  5  . 

a  function 

,  as 

shown 

in 

Figure 

E(x)  =  max 


(x  -  b),  (b  -  x) 


9  o 


»  o 


=  max  (u+  e+,  u'  e*)  . 

Notice  that  the  penalty  value  will  oe  relatively  small  (less  than  1 
values  between  b_  and  F  and  that  it  will  increase  at  a  rate  1  / ( b  ■*) 
x  <  b  or  l/(b  9)  if  b  <  x. 


1  for 

if 
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E(x) 


Fig.  1.  Piecewise  Function  for  Linear  Programming. 


3,  Elastic  Programming 

For  this  type  of  programming  model,  the  definition  of  the  error 
will  be  modified  slightly  as  shown  in  Taole  7.  The  weighting  factors 
will  be  changed  since  the  actual  error  is  now  slightly  different  due  to 
the  use  of  the  error  bounds  rather  than  the  actual  requirements.  In 
order  to  notationally  indicate  this  difference,  the  use  of  underscore  and 
overscore  for  weights  and  errors  will  be  used  in  place  of  superscript 
plus  and  minus.  This  type  of  programming  can  use  a  piecewise  linear 
objective  function  of  the  following  type: 

.(b  -  x) 


E(x)  =- 


*  b 
0 

(x  -  F) 
9  F 


if  x  <  b_  , 
if  b_  £  x  <_  F  , 
if  F  <  x  . 


=  max  (u  e,  0,  u  e)  . 

This  function  is  shown  in  Figure  2.  Notice  that  the  penalty  is  0  for 
values  of  x  within  the  range  of  permitted  errors  (_b,  o’)  and  increases  at 
a  rate  l/{o  *)  if  x  is  less  than  D_  or  1/(F  9)  if  x  is  greater  than  F.  Also, 
notice  that  if  b  =  b  =  F  then  this  objective  is  the  same  as  the  goal 
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TABLE  7 


ERROR  FOR  ELASTIC  PROGRAMMING 

by  job  type  by  people  type 

and  rank  and  rank 


shortf al 1 

lik 

'  -2-ik 

*  Zik 

overf i 1 1 

eu  ■  xij  -  blj 

*ik 

N 

II 

*  °ik 

E(x) 

x 

b.  b 

Fig.  2.  Piecewise  Function  for  Elastic  Programming 

programming  objective.  Thus,  the  goal  program  is  a  special  case  of  the 
elastic  program. 

4.  Quadratic  °roqramminq 

Several  quadratic  programming  formulations  are  presented  in 
Section  i/I.  These  formulations  are  easy  to  use  and  give  quick  results 
for  formulations  that  can  oe  expressed  without  inequality  constraints. 

The  quadratic  objective  function  to  be  used  must  have  the  same  properties 
mentioned  *'or  the  piecewise  ’inear  function  and  it  must  also  be  greater 
than  zero.  These  properties  may  be  stated  mathematically  as: 


(i) 

E(x) 

=  0 

if  x  =  b  , 

lii) 

E  ( x ) 

=  1 

if  x  =  b  , 

(iii) 

E(x) 

=  1 

if  x  =  F  . 

( iv) 

E(x) 

>  0 

> 

Unfortunately,  quadratic  functions  which  use  the  previously  defined 
penalty  measure  must  have  symmetric  penalties  in  order  to  satisfy  the  four 
conditions.  That  is,  y  must  equal  9  and  5  must  equal  *  when  considering 
the  decision  variable  z.  The  following  objective  satisfies  these 
conditions: 

e(x)  =  (2Lir)2  when  9  =  *  • 

This  type  of  penalty  function  is  shown  in  Figure  3.  An  approximation 
that  is  useful  if  9  and  y  are  close  but  not  equal  is  described  in 
Grinold  [Ref.  1]. 


E(x) 

I 

1 

I 

1 

I 

1  | 


b  b  b 

Fig.  3.  Function  for  Quadratic  Programming. 

£ .  PROGRAMMING  NOMENCLATURE 

It  will  be  useful  to  have  a  notational  scheme  to  identify  the 
various  orogramming  types,  decision  variables,  and  objective  functions 
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found  in  the  following  sections.  Let  LP,  EP,  and  QP  refer  to  linear 
(goal),  elastic,  and  quadratic  programming,  respectively.  Let  the 
objective  function  that  sums  errors  by  job  class  be  a  type  1  objective 
function,  the  objective  function  that  sums  the  errors  by  people  class  be 
type  2,  and  the  one  that  sums  errors  by  both  job  class  and  people  class 
be  type  3.  Thus,  the  mathematical  program  LP1(Y)  is  a  linear  program 
that  uses  the  first  objective  function  with  Y  as  a  decision  variable. 
Similarly,  QP2(Y,X)  is  a  quadratic  program  with  the  second  type  of 
objective  function  that  uses  Y  and  X  as  the  decision  variables. 


F.  ADDING  COST  CONSIDERATIONS 

Grinold  [Ref.  1]  also  discusses  a  way  in  which  cost  considerations 
may  be  added  to  the  model.  Let  c^  be  the  annualized  cost  of  a  type  k 
accession.  Then  the  multi-attribute  objective  function: 


X  E  C 
k 


k  v+  ^  \z  z  Eii  +  2 

i  j  J  i 


E  £  0: 


ik ' 


where  0  £  x  <_  1,  will  minimize  the  objective  with  respect  to  both  cost 
and  weighted  penalties.  This  objective  leads  to  a  family  of  objective 
functions  depending  on  the  value  of  the  parameter  x.  Notice,  however, 
that  the  units  of  the  objective  function  are  now  in  terms  of  cost  units 
and  percent  error  units. 

The  use  of  such  multi-attribute  functions  is  advantageous  when 
the  annualized  cost  of  one  type  of  manpower  is  significantly  different 
than  another  type.  For  example,  the  annualized  cost  of  a  general  unre¬ 
stricted  line  officer  is  vastly  less  than  that  of  a  pilot  aue  to  traininq 
costs,  oay  bonuses,  etc. 
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IV.  LINEAR  (GOAL)  PROGRAMMING  FORMULATION 


A.  GENERAL 

A  linear  programming  formulation  for  minimization  of  errors  is 
found  in  Grinola  [Ref.  1]  and  is  derived  from  the  goal  programming 
concept  as  described  in  Lee  [Ref.  5]  and  Hillier  and  Lieberman  [Ref.  6]. 
Goal  programming  is  a  modification  and  extension  of  linear  programming 
that  allows  a  simultaneous  solution  of  a  system  of  complex  objectives 
rather  than  a  single  objective.  Notice  that  the  goal  of  minimizing 
errors  by  differing  ranks,  job  types,  and  people  types  constitutes  a 
system  of  conflicting  objectives.  A  solution  which  forces  exact  achieve¬ 
ment  of  goals  for  the  lower  ranks  may  result  in  extreme  deficits  in  the 
higher  ranks  while  one  which  exactly  meets  requirements  in  the  higher 
ranks  may  result  in  large  excesses  in  the  lower  ranks.  Thus,  the  goal 
approach  is  used  to  seek  the  best  possible  solution  given  the  stated 
trade-offs. 

In  implementation  of  the  following  programs,  the  user  should  be  aware 
of  the  tactical  problems  in  solving  LP's.  It  may  be  advantageous  in 
some  cases  to  solve  the  dual  formulations  of  the  problems  depending  on 
the  number  of  variations  and  the  number  of  constraints.  An  example  of 
forming  the  dual  for  LP3( Y }  and  LP3(X,Y)  is  in  Grinold  [Ref.  1]  and  a 
general  discussion  of  the  topic  is  in  Luenberger  [Ref.  7], 

3.  ACCESSION  ANO  INVENTORIES 

The  first  formulation  will  be  designated  LP1(Y)  and  is  shown  in 
Figure  4.  Recall  that  lPL  formulations  involve  errors  of  rank  and 


job  type.  The  constraint  arises  from  the  conservation  of  flow  derived 
from  equation  (6)  and  sharing  property  in  equation  (7).  An  alternate 
formulation  that  considers  both  types  of  errors  is  LP3(Y)  as  shown  in 
Figure  5.  This  is  labelled  in  Grinold  [Ref.  1]  as  LP-P. 
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Figure  6  displays  the  formulation  which  uses  both  X  and  Y  as  the  decision 
variables.  It  is  LP3(X,Y)  and  is  round  in  Grinold  [Ref/  1]  as  LP-B. 


C.  EXPECTED  STAGE  LENGTHS 

Another  application  using  qoal  programming  would  be  to  minimize 
the  errors  with  respect  to  expected  stage  lengths  (W).  Recall  that  the 
stage  length  for  rank  i,  person  type  k  is  the  expected  amount  of  time  a 
recruit  in  class  k  will  spend  in  rank  i.  An  added  constraint  for  this 
formulation  is  that  the  stage  lengths  have  upper  and  lower  bounds.  In 
the  Navy  officer  example,  such  bounds  would  reflect  navy  policy  or  legal 
restrictions  imposed  by  Congress.  Let  these  upper  and  lower  bounds  be 
defined  by  matrices  W  and  W,  respectively.  LP1(W)  may  be  expressed  as 
shown  in  Figure  7. 
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Fig.  7.  LPl(W) . 


0.  SHARING  RULES 

Of  particular  interest  is  a  mathematical  program  to  specify  the 
sharing  fraction  that  represents  the  policy  for  distributing  job  classes 
among  the  people  types  (F)  or  people  classes  among  job  types  (G).  An 
additional  constraint  on  this  formulation  is  that  F  must  sum  to  one  over 
the  index  k  and  G  must  sum  to  one  over  the  index  j.  Also,  the  sharing 
fractions  may  be  bounded  above  and  below.  Using  the  same  convention  as 
used  in  the  case  of  stage  lengths,  let  F_  and  _G  be  the  lower  bounds  and 
let  T  and  GT  be  the  upper  bounds  on  F  and  G,  r.  pectively.  Notice  that  if 
no  sharing  is  possible,  specifying  that  both  the  upper  and  lower  bound  be 
zero  sets  the  value  at  zero.  LP1(G)  and  LP2(F)  are  formulated  as  shown 
in  Figures  3  and  9. 
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Fig.  8.  LP1(G) 
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E.  INTERPRETATION  FROM  LP  THEORY 

In  classical  LP  terms,  the  errors  e+,  e* ,  d+,  and  d‘  represent  logical 
variables  introduced  to  maintain  the  equality  constraints.  The  errors  e* 
and  d"  may  be  thought  of  as  artificial  variables  and  the  errors  e+  and  d+ 
as  surplus  variables  and  hence,  they  measure  the  distance  from  a  particular 
solution  to  the  requirement.  These  objective  functions  are  mathematically 
equivalent  to  a  Lagrangian  form  with  cost  coefficients  of  the  decision 
variables  equal  to  zero  (see  Duff  [Ref.  8:  p.  64]).  The  weighting  factors 
u+,  u",  v+,  and  v"  represent  bounds  on  the  variables  of  the  dual  formulation. 
In  fact,  an  interpretation  of  a  particular  weighting  factor  is  that  it 
defines  the  penalty  cost  per  unit  violation  of  its  associated  constraint. 
Therefore,  weighting  factors  which  are  approximately  equal  have  associated 
constraints  which  are  of  the  same  relative  importance,  if  binding.  In 
other  words,  the  incremental  penalty  for  violating  each  such  constraint 
is  approximately  the  same.  Such  information  should  be  of  use  in  analyzing 
the  results  of  these  programs. 
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V.  ELASTIC  PROGRAMMING  FORMULATION 


A.  GENERAL 

The  elastic  programming  model  is  described  in  Duff  [Ref.  8].  The 
formulation  is  quite  similar  to  that  of  the  goal  program  in  the  last 
section.  However,  there  is  now  a  subtle  but  important  difference.  In 
the  goal  program,  the  objective  was  to  allocate  people  as  closely  as 
possible  to  an  exact  number  of  billets  (either  8  or  P).  In  the  elastic 
program,  the  objective  is  now  to  allocate  people  so  as  to  be  within  a 
specified  range  of  the  billets  to  be  filled.  For  this  reason,  the  actual 
error  will  be  indicated  by  a  change  in  the  notation  to  the  use  of  under¬ 
score  and  overscore  in  place  of  superscripts  plus  and  minus. 

Before  starting  the  formulations ,  it  is  necessary  to  define  a  variable 
to  represent  the  slack  from  the  upper  bounds  (3  or  P).  Let  and  r 
represent  these  slack  variables  for  the  (i,j)  and  (i,k)  inventories.  Also 
recall  that  _u,  IT,  v_,  v  are  the  weighting  factors  and  e_,  e,  _d,  d-  are  the 
errors  from  the  desired  ranges. 

8.  ACCESSION  AND  INVENTORIES 

The  first  formulation  shown  in  Figure  10  below  is  similar  to  the 
first  formulation  in  the  previous  chapter  and  is  denoted  EP1(Y).  The 
variable  Y  once  again  represents  an  upper  bound  on  the  accession  vector  Y. 
An  alternate  formulation  that  considers  both  type  of  errors  is  EP3(Y) 
presented  in  Figure  11. 
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C.  EXPECTED  STAGE  LENGTHS  AND  SHARING  RULES 


Elastic  programming  may  also  be  used  to  Determine  expected  stage 
lengths  or  sharing  rules.  Only  one  example,  EP2 ( F )  in  Figure  12,  will  be 
presented  since  the  formulations  are  quite  similar  to  those  in  the  last 
section.  Recall  that  an ’additional  constraint  exists  for  the  sharing 
rules  since  they  must  sum  to  one  over  the  appropriate  index.  An  example 
of  this  is  £P2(F)  shown  below. 


Fig.  12.  EP2(F) 

D.  INTERPRETATION  FROM  LP  THEORY 

As  was  the  case  for  the  linear  (goal)  program,  the  elastic  program 
is  mathematically  equivalent  to  a  Lagrangian  form.  The  variables  £  and  d 
are  artificial  and  measure  the  distance  from  a  particular  solution  to  the 
lower  error  bound.  The  surplus  variables  ?  and  <T  measure  the  distance 
above  the  upper  error  bound  and  the  slack  variables  q  and  r  measure  the 
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distance  from  a  solution  to  the  associated  upper  bound.  The  weighting 
factors  continue  to  represent  bounds  on  the  dual  formulation. 

This  formulation  has  an  advantage  in  that  the  solution  space  for 
zero  penalty  is  greatly  expanded  even  if  the  error  range  is  small.  This 
has  improtant  consequences  for  the  implementation  and  achievement  of 
acceptable  solutions  for  problems  with  many  constraints  and  variables. 
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VI.  QUADRATIC  PROGRAMMING  FORMULATIONS 


A.  GENERAL 

The  quadratic  programming  model  described  in  Grinold  [Ref.  1]  pro¬ 
vides  an  alternate  method  for  optimization  of  allocations.  A  general 
discussion  of  the  solution  techniques  for  quadratic  programs  is  found  in 
Simmons  [Ref.  9].  The  following  programs  are  expressed  with  only  equality 
constraints  and  will  yield  solutions  that  are  found  by  analytically 
solving  systems  of  equations.  This  provides  a  useful  and  consistent 
means  of  determining  allocations.  However,  the  use  of  QP  as  presented 
here  must  be  weighed  against  two  disadvantages. 

•  The  penalty  weights  must  De  symmetric  for  both  deficit  and  excess 
types  of  errors. 

•  The  decision  variables  are  unbounded  and,  thus,  the  QP  may  yield 
impractical  solutions  such  as  negative  accessions. 

Since  much  of  the  data  input  for  the  model  is  subjective  or  highly 
variable,  and  the  purpose  of  the  model  is  to  explore  questions  concerning 
policy  planning,  the  magnitudes  and  relative  relationships  of  the  alloca¬ 
tion  are  more  important  than  an  exact  choice  of  allocation.  Thus,  the 
use  of  symmetric  penalties  may  not  present  significant  problems.  If  the 
parameters  and  data  inputs  to  the  model  are  fairly  consistent  with  real 
world  conditions,  then  the  use  of  unbounded  variables  is  unlikely  to 
produce  inconsistent  solutions  which  would  adversely  affect  tne  result  of 
the  analysis.  The  appearance  of  negative  values  in  the  solution  Is 
likely  to  indicate  misspecif ication  of  parameters  or  the  attempt  to 
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achieve  unrealistic  goals.  For  the  following  sections,  it  will  be 
assumed  that  ^  =  0  and  s  =  <*>  . 


The  QP  is  similar  to  a  least  squares  approach  in  that  the  purpose  is 
to  minimize  the  sum  of  squared  weighted  errors.  Appendix  B  contains  a 
brief  summary  of  the  resulting  quadratic  forms  used  in  this  section. 
Reindexing  the  variables  makes  it  possible  to  conveniently  express  the  QP 
as  a  system  of  equations  using  matrix  arithmetic  notation.  For  mathemat¬ 
ical  convenience,  the  following  change  of  variables  may  be  necessary: 

m  =  J  (i  -  1)  +  j  (8) 

where  m=l,2,...M  and  M=IJ; 

n  ■  K  (i  -  1)  +  k  (9) 

where  n«l,2,...,N  and  N=IK. 

In  this  way,  matrices  of  size  (I,K)  may  be  reindexed  to  become  vectors 
of  size  (N)  and  matrices  of  size  (I,J)  become  vectors  of  size  M.  Three 
dimensional  arrays  may  also  be  reduced  to  two  dimensions.  For  example,  a 
three  dimensional  array  of  size  (I,J,K)  indexed  over  (i,j,k)  may  be 
reindexed  as  a  two  dimensional  array  of  size  (M,K)  using  equation  (8). 

In  using  matrix  notation,  it  will  be  necessary  on  occasion  to  reform  a 
vector  into  a  diagonal  matrix.  This  means  that  the  kth  element  of  the 
vector  will  occupy  the  (k,k)th  position  of  the  matrix  and  all  off 
diagonal  elements  will  be  zero. 

3.  ACCESSION  AND  INVENTORIES 

The  first  formulation  to  be  considered  is  QP1(V).  The  advantage  of 
this  formulation  is  that  it  requires  less  data  input  since  the  job 
sharing  fractions  (F)  and  requirements  by  rank  and  manpower  class  (P)  3re 
not  needed.  Also,  once  the  optimal  accession  and  inventories  are 
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determined,  values  for  F  and  P  may  be  computed  based  on  the  optimization 


results.  The  objective  function  is: 


min  z  Z 
i  J 


uij(  k2  gkij  *ik  \  '  bij  > 


ij 


Now  let  h^  =  gk  •  j  wik  and  then  reindex  (i,j)  to  (m)  and  (i,j,k)  to 
(m,k)  as  given  above  in  equation  (8).  This  results  in: 


u2 

m 


min  £  u1,  (  2  h  ,  v,  -  b  ) 
-  '  mk  -  k  nr 


m  k 

which,  in  matrix  notation,  equals: 


min  (HY  -  B)1  U’U  (HY  -  B) 

where  U  is  now  an  (M,M)  diagonal  matrix  of  the  penalty  weights.  Taking 
the  derivative  with  respect  to  Y  and  setting  the  result  equal  to  zero 
yields  the  following  solution: 

Y  =  (H'U'UH)*1  H'U'UB  . 


Once  an  optimal  value  for  Y  has  been  determined,  inventories  may  be 
calculated  using  equations  (7)  and  (6)  from  Section  II. 

The  formulation  QP3(Y)  minimizes  the  objective: 

*in  f  j  4(  xu  -  v2  *  f  k*  4  <"ik  *k  •  bik>2 

suoject  to  the  sharing  constraint 


Xi j  *  k  gk i j  wik  yk  * 

Substituting  the  constraint  into  the  objective,  letting  h.jk  =  gki^  w-k, 
ana  then  reindexing  from  (i,j,k)  to  (m,k)  as  given  above  in  equation  (8), 
yields  in  matrix  notation: 

min  (HY  -B)‘  U'U  (HY  -  3)  +  Y’QY  -  2Y'R  +  c 
where  Q  is  a  (K,K)  diagonal  matrix  whose  diagonal  elements  are  given  by 
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2 

R  is  a  (K)  vector  with  elements  rk  =  z  v^k  p-k  wik  ; 

2  2 

and  c  is  a  scalar  constant,  namely  c  =  i  s  v.,  p..  . 

i  k  1K 

Taking  the  derivative  and  setting  it  equal  to  zero  results  in 
H'U'UHY  -  H'U'UB  +  QY  -  R  =  0  . 

Finally,  the  solution  is: 

Y  =  [H'U'UH  +  Q]'1  [H'U'UB  +  R]  . 

This  model  can  be  found  in  Grinold  [Ref.  1]  as  UQ-P.  An  advantage  to 
the  use  of  this  model  over  other  QP3  formulations  is  that  it  uses  the 
people  sharing  rule  G  which  may  be  more  available  as  a  data  input  than 
the  job  sharing  rule  F. 

An  alternative  program  which  uses  the  same  objective  function  is 
QP3( X, Y) ,  In  this  case,  minimize  the  objective  subject  to  the  sharing 
constraint: 


"ik  -  s  fkij  xij  *  0  for  3,1  '•k- 

J 

Both  X  and  Y  are  now  the  decision  variables  and  the  constraint  cannot 
be  substituted  into  the  objective  function  directly.  In  order  to  exDress 
the  constraint  in  matrix  notation,  it  is  necessary  to  define  two  new 
matrices,  F*  and  W*.  First,  F*  is  defined  by  the  elements: 


fi '  ki  j 


f.  .  .  if  i '  =  i 
kij 


0  otherwise 
where  i,  i'=l,...,I,  k=l,...,K,  and  j=I,...,J. 

Now  reindex  from  (i',k,i,j)  to  (n,m)  as  given  by  equations  (8)  and  (9) 
above  to  form  the  properly  shaped  F*.  This  results  in  a  block  matrix  in 
which  each  of  the  blocks  is  a  diagonal  matrix  of  the  jth  column  of  F 
for  a  particular  (i,k).  W*  is  formed  in  a  similar  manner.  Let: 
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for  k  =  k' 


1^0  otherwise 
where  k,  k'=l,...,K  and 

Next  reindex  W*  from  (i,k',k)  to  (n,k)  according  to  (9)  above. 

The  result  is  a  diagonal  block  matrix  in  which  each  of  the  blocks  on 
the  diagonal  is  a  column  of  W.  Finally,  the  constraint  can  be  expressed 
in  matrix  notation  as: 

F*X  -  W*Y  =  0  . 

The  objective  function  in  matrix  notation  is: 

min  (X  -  3)'  U'U  (X  -  3)  +  Y'QY  -  2Y'R  +  c 
where  1),  Q,  R,  and  c  are  the  same  as  defined  for  QP3(Y).  Now  form  the 


Lagrangian: 


(X  -  B) 1  U'U  (X  -  B)  +  Y'QY  -  2Y'R  +  c  -  x  (F*X  -  W*Y) 


where  x  is  the  vector  of  Lagrange  multipliers.  Taking  the  derivative 
with  respect  to  X,  Y,  and  x  yields  the  following  system  of  equations: 
2U'UX  -  2U'UB  -  F*'  ■  0, 

2QY  -  2R  +  W*  =  0, 

F*X  -  W*Y  =  0. 

This  can  be  written  in  partitioned  matrix  form  as: 
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where  0  indicates  a  zero  matrix  of  appropriate  size.  The  optimal  Y,  X, 


and  \  is  given  by: 
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The  program  QP3 ( X , Y ) 

is  found 

i  n 

Grinold  [Ref.  1]  with 

the 

name  UQ 

can  be  seen  from  the  result,  QP3(X,Y)  has  a  more  complicated  solution 
than  QP3(Y)  and  has  an  additional  disadvantage  in  that  the  results  may 
slightly  violate  the  consistency  rule.  Thus,  a  determination  must  be 
made  whether  to  use  both  X  and  Y  from  the  results  or  to  construct  the 
inventories  solely  on  the  basis  of  Y.  It  has  been  presented  here  as  an 
illustration  of  the  use  of  Lagrange  multipliers  and  as  a  means  of  checking 
the  results  from  QP3(Y). 

Another  means  of  checking  QP3 ( Y )  which  insures  that  the  consistency 
rule  is  met  would  be  to  solve  QP3(X),  use  equation  (5)  to  determine  z, 
and  then  use  the  least  squares  approach  presented  in  Section  II  to  find 
values  for  Y.  Those  values  of  Y  may  then  be  used  to  reconstruct  the 
inventories  and  insure  consistency. 

C.  EXPECTED  STAGE  LENGTHS 

The  next  formulation  that  will  be  discussed  is  QP2(W).  The  objective 


function  is: 


”1n  1  ;  vik  (»k  “ik  •  pik>2 

w  t  k 

Let  qik  =  vik  *k  and  rik  =  vik  pik  ' 

Then  the  result  is: 

min  r  £  ( q i k  wik  *  rik)2‘ 
w  i  k 

Now  reindexing  from  (i,k)  to  (n)  and  changing  to  matrix  notation  yields: 
min  (QW  -  R) '  (QW  -  R) 

where  Q  is  reformed  to  be  an  (n,n)  diagonal  matrix.  The  solution  to 
this  program  is: 

W  =  [Q'Q]"1  Q'R 

D.  ADDING  ADDITIONAL  EQUALITY  CONSTRAINTS 

The  technique  of  Lagrange  multipliers  presented  with  the  QP3(X,Y) 
model  may  also  be  used  with  additional  equality  constraints.  For  example, 
suppose  the  total  number  of  accessions  to  be  brought  into  the  system  is 
Known  but  the  exact  number  in  each  manpower  class  is  not  known.  This 
information  may  be  added  as  an  additional  equality  constraint: 

£  Yu  =  r 
k  K 

where  r  is  now  the  total  number  of  recruits  to  be  brought  into  the 
system  each  time  period.  In  matrix  notation,  this  constraint  becomes : 
e'Y  -  r  =  0 

where  e  is  a  vector  of  k  ones.  Using  the  QP1(Y)  model  as  an  example, 
form  the  Lagrangian  with  this  additional  constraint  to  get 
min  (HY  -  3)'  !J'U  (HY  -  3)  -  X  ^  e '  v  -  r)  . 
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Take  the  derivatives,  set  equal  to  zero,  and  solve  as  before  to  get 
the  solution: 


r 

1 

1 

r 

1 

T1 

; 

r  -i 

1  Y 

1 

1 

| 

|  H'U'UH 

! 

|  | 

1  -{1/2)1  1 

1  1 

|  | 

1  H'U'UB  | 

1  1 

i  _  i 

|  X 

1 

1 

I 

1 

1  e  1 

i 

- -  j 

!  1 

i  -r  | 

I  | 

1  i  : 

10  1 

i  i  i 

L 

I 

J. 

1 

1  1 

1  J 

r 

i  _ 

ii  I.  .'rfiii 

The  advantage  of  adding  an  additional  constraint  is  that  it  can  be 
used  to  force  the  quadratic  model  to  search  for  solutions  that  will  be 
more  related  to  "real  world"  problems.  This  helps  to  mitigate  the 
disadvantage  regarding  unbounded  decision  variables.  In  fact,  a  range  of 
possible  values  could  be  specified  and  then  the  optimal  solution  within 
the  range  could  be  found  using  a  simple  numerical  search  procedure. 

Since  the  solution  to  each  problem  during  a  search  iteration  is  explicit 
and  the  resolution  needed  is  to  the  nearest  whole  integer,  such  a  proce¬ 
dure  is  easy  to  implement  and  should  arrive  at  solutions  fairly  quickly. 
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VII.  DEMONSTRATION  OF  THE  MODEL 


A.  IMPLEMENTATION 

The  model  was  implemented  in  its  quadratic  form  on  an  IBM  3033 
computer  at  the  Naval  Postgraduate  School  using  the  APL  programming 
language  as  shown  in  Appendix  C.  For  the  linear  and  elastic  programs, 
preprocessors  were  written  in  APL  to  arrange  the  data  in  standard  MPS 
programming  format  for  implementation  on  the  XS-system  optimization 
package.  MPS  is  an  international  standard  and  such  format  is  easily 
transferrable  to  other  types  of  computers.  A  sample  output  is  included 
in  the  program  listing  section.  In  addition,  listings  of  some  of  the  APL 
program  functions  used  to  produce  the  results  in  this  section  are  also 
included.  The  ease  with  which  data  arrays  can  be  manipulated  in  APL 
resulted  in  concise  and  flexible  program  functions  which  performed  the 
optimizations  and  the  input-output  formatting  of  the  data.  Also  written 
were  two  "user  friendly"  functions  for  data  input  and  data  display,  and  a 
master  driver  function  that  allows  the  user  to  select  from  a  menu  of 
options.  Appendix  C  contains  a  brief  summary  of  these  programs. 

8.  SYSTEM  DESCRIPTION 

1.  Attributes 

This  demonstration  is  based  on  an  example  from  the  U. S.  Navy 
Officer  Corps.  The  attributes  chosen  for  manpower  classes,  career  stages, 
ana  job  classes  are  the  same  as  those  found  in  Grinold  [Ref.  1]  3nd  are 
shown  in  Figure  13. 
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MANPOWER  CLASSES 


Type  Officer  Type 


Description 


1  llOx 

2  11 lx / 1 16x 

3  112x/117x 

4  13  lx / 1 39x 

5  132x/137x 


General  Unrestricted  Line  (GURL) 
Surface  Warfare 
Submarine  Warfare 
Pilots 

Naval  Flight  Officer  (NFO) 


CAREER  STAGE 


Stage 


Rank  Years  of  Service 


1  ENS  0-2 

2  LTJG  2  -  4 

3  LT  4-9 

4  LCDR  9  -  14 

5  CDR  14  -  19 

6  CAPT  19  -  26 


JOB  CLASSES 


Type  Billet  Designator 


Description 


1  1000 

2  1050 

3  1110/1160 

4  1120/1170 

5  1310/1390 

5  1320/1370 

7  1300 


General,  Nonwarfare 
General,  Warfare 
Surface  Warfare 
Subsurface  Warfare 
Pilots 

Naval  FI ight  Officer 
General  Aviation 


Fig.  13.  Attributes  for  Example. 


2.  Expected  Time  in  Each  Rank 

To  compute  W,  it  is  necessary  to  have  stage  lengths  and  survivor 
fractions.  Navy  usage  of  survivor  fractions  is  usually  found  in  the  form 
of  continuation  or  retention  rates.  Current  data  was  obtained  +'rom 
Deputy  Chief  of  Naval  Operations  (MPT),  OP -01,  with  continuation  '-ates 
for  Na\y  Unrestricted  Line.  The  continuation  rates  were  used  directly  to 
compute  expected  length  of  time  spent  in  each  rank.  This  data  is  snown  in 
Appendix  D.  The  APL  function  WAITS  uses  continuation  rates  and  staae 


53 


lengths  as  inputs  and  produces  the  matrix  W.  the  program  GWAITS  will 
also  accomplish  this  but  uses  survivor  fractions  instead  of  retention 
ra  .  The  resulting  steaay  state  matrix  W  is  shown  in  Table  8. 


TABLE  8 

EXPECTED  STAGE  LENGTHS 


GURL 

SURF 

SUB 

PILOT 

NFO 

ENS 

1.925 

1.942 

1.891 

2.000 

1.992 

LTJG 

1.674 

1.700 

1.621 

1.975 

1.887 

LT 

2.063 

2.073 

2.040 

2.827 

3.241 

LCDR 

1.106 

1.138 

1.040 

1.028 

2.255 

CDR 

0.856 

0.896 

0.777 

0.773 

1.869 

CAPT 

0.596 

0.589 

0.607 

0.445 

1.276 

3.  Billet  Requirements  and  Sharing  Rules 

Information  concerning  the  requirements  by  rank  and  job  type  was 
also  obtained  from  0P-01  and  the  array  B  was  constructed.  (See  Table 
9).  The  next  parameter  to  be  specified  is  the  billet  sharing  array.  This 
information  was  taken  from  Grinold  [Ref.  a]  with  some  slight  modifications 
and  is  displayed  in  Table  10.  Once  B  and  F  are  available,  it  is  possible 


TABLE  9 


BILLET  REQUIREMENTS  BY  RANK  AND  JOB  TYPE 


1000 

1050 

1110 

1120 

1310 

1320 

1300 

ENS 

330 

0 

2258 

660 

1105 

593 

0 

LTJG 

768 

0 

1706 

700 

1970 

1185 

2 

LT 

1780 

378 

2145 

884 

3780 

1408 

629 

LCDR 

1456 

571 

1510 

377 

1818 

586 

716 

CDR 

1001 

468 

990 

518 

710 

54 

909 

CAPT 

571 

462 

354 

152 

0 

J 

422 

54 


- —  - 


rdrr'ni-rffcii 


TABLE  10 


BILLET  SHARING  ARRAY 


GURL 

1000 

1050 

1110 

1120 

1310 

1320 

1300 

ENS 

0.73 

0 

0 

0 

0 

0 

0 

LTJG 

0.70 

0 

0 

0 

0 

0 

0 

LT 

0.70 

0 

0 

0 

0 

0 

0 

LCDR 

0.70 

0 

0 

G 

0 

0 

0 

CDR 

0.70 

0 

0 

0 

0 

0 

0 

CAPT 

0.70 

0 

0 

0 

0 

0 

0 

SURF 

1000 

1050 

1110 

1120 

1310 

1320 

1300 

ENS 

0.18 

1.00 

1 

0 

0 

0 

0 

LTJG 

0.15 

0.50 

1 

0 

0 

0 

0 

LT 

0.15 

0.50 

1 

0 

0 

0 

0 

LCDR 

0.15 

0.50 

1 

0 

0 

0 

0 

CDR 

0.15 

0.50 

1 

0 

0 

0 

0 

CAPT 

0.15 

0.50 

1 

0 

0 

0 

0 

SUB 

1000 

1050 

1110 

1120 

1310 

1320 

1300 

ENS 

0.09 

0.00 

0 

1 

0 

0 

0 

LTJG 

0.06 

0.10 

0 

1 

0 

0 

0 

LT 

0.06 

0.20 

0 

1 

0 

0 

0 

LCDR 

0.06 

0.20 

0 

1 

0 

0 

0 

CDR 

0.06 

0.20 

0 

1 

0 

0 

0 

CAPT 

0.06 

0.20 

0 

1 

0 

0 

0 

PILOT 

1000 

1050 

1110 

1120 

1310 

1320 

1300 

ENS 

0.00 

0.00 

0 

0 

1 

0 

0.00 

LTJG 

0.50 

0.10 

0 

0 

1 

0 

0.00 

LT 

0.50 

0.18 

0 

0 

1 

0 

0.57 

LCDR 

0.50 

0.18 

0 

0 

1 

0 

0.57 

CDR 

0.50 

0.18 

0 

0 

1 

0 

0.57 

CAPT 

0.50 

0.18 

0 

0 

i 

X 

0 

0.57 

NFO 

1000 

1050 

1110 

1120 

1310 

1320 

1300 

ENS 

0.00 

0.00 

0 

0 

0 

1 

1.00 

LTJG 

0.40 

0.10 

0 

0 

0 

1 

1.00 

LT 

0.40 

0.12 

0 

0 

0 

1 

0.43 

LCDR 

0.40 

0.12 

0 

0 

0 

1 

0.43 

COR 

0.40 

0.12 

0 

0 

0 

1 

0.43 

CAPT 

0.40 

0.12 

0 

0 

0 

1 

0.43 
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to  calculate  the  corresponding  requirements  in  terms  of  people  types  and 
ranks,  and  the  corresponding  people-sharing  array.  The  results  for  the 
requirements  by  people  types  is  shown  in  Table  11.  An  example  of  the 
people-sharing  array  is  shown  in  Table  12  for  the  manpower  class  of 
pilots. 


TABLE  11 

BILLET  REQUIREMENTS  BY  RANK  AND  MANPOWER  CLASS 


GURL 

SURF 

SUB 

PILOT 

NFO 

ENS 

241 

2317 

690 

1105 

593 

LTJG 

538 

1821 

746 

2010 

1217 

LT 

1246 

2601 

1066 

4296 

1795 

LCDR 

1019 

2014 

1079 

2402 

1021 

CDR 

701 

1374 

672 

1362 

541 

CAPT 

400 

671 

279 

352 

260 

TABLE  12 

PEOPLE-SHARING  ARRAY  FOR  PILOTS 


1000 

1050 

1110 

1120 

1310 

1320 

1300 

ENS 

0.015 

0.000 

0 

0 

0.985 

0 

0.000 

LTJG 

0.019 

0.000 

0 

0 

0.980 

0 

0.001 

LT 

0.021 

0.016 

0 

0 

0.880 

0 

0.083 

LCDR 

0.030 

0.043 

0 

0 

0.757 

0 

0.170 

CDR 

0.037 

0.062 

0 

0 

0.521 

0 

0.380 

CAPT 

0.081 

0.236 

0 

0 

0.000 

0 

0.683 

4.  Permitted  Unit  Errors 

The  permitted  unit  errors  were  determined  based  on  subjective 
judgment.  First  the  matrix  for  permitted  unit  errors  by  rank  and  job 
type  was  estimated.  This  is  shown  in  Table  13.  The  entries  marked  as  1 
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TABLE  13 


PERMITTED  UNIT  ERRORS  BY  JOB  TYPE  AND  RANK 


1000 

1050 

1110 

1120 

1310 

1320 

1300 

ENS 

400 

1 

500 

500 

500 

400 

1 

LTJG 

400 

1 

500 

500 

500 

400 

500 

LT 

10 

3 

8 

10 

8 

10 

10 

LCDR 

10 

3 

10 

15 

3 

8 

10 

CDR 

10 

8 

5 

10 

8 

10 

10 

CAPT 

10 

10 

10 

10 

1 

1 

10 

correspond  to  ranks  and  job  types  for  which  there  are  no  billets. 

Giving  a  very  low  percentage  will  help  ensure  the  optimization  does  net 
seek  to  place  people  in  these  jobs.  The  corresponding  errors  by  people 
type  and  rank  were  generated  by  the  use  of  the  consistency  rule.  The 
unit  errors  do  not  necessarily  have  to  meet  this  requirement,  but  it  was 
felt  that  the  consistency  rule  would  give  an  approximation  as  reasonable 
as  further  use  of  subjective  judgment.  The  results  were  checked  to  insure 
that  the  numbers  were  reasonable,  and  no  changes  were  made.  The  matrix 
for  permitted  unit  errors  by  rank  and  people  type  is  shown  in  Table  14. 
These  last  two  matrices  represent  the  permitted  error  for  overfilling  the 
billets.  The  matrices  for  underfilling  the  billets  were  done  using  a 
simi 1 ar  method . 

C.  ACCESSIONS  USING  QUADRATIC  PROGRAMMING 

The  first  set  of  programs  to  be  considered  are  those  labelled  QP1(Y), 
QP2(Y),  and  QP3(X,Y).  The  results  of  the  optimization  with  respect  to 
the  accessions  are  shown  in  Table  15.  The  column  marked  "Obj.  Val."  is 
the  value  of  the  penalty  function  that  measures  both  types  of  errors. 
Notice  that  QP1(Y)  wnich  only  minimized  with  respect  to  errors  by  job 


TABLE  14 


PERMITTED  UNIT  ERRORS  BY  RANK  AND  PEOPLE  TYPE 


GURL 

SURF 

SUB 

PILOT 

NFO 

ENS 

280.0 

560.5 

524.2 

520.7 

416.5 

LTJG 

280.0 

560.5 

524.2 

805.2 

631.1 

LT 

7.0 

13.5 

12.2 

15.6 

15.7 

LCDR 

7.0 

15.5 

17.2 

15.6 

13.7 

CDR 

7.0 

10.5 

12.2 

15.6 

15.7 

CAPT 

7.0 

16.5 

12.6 

9.0 

6.9 

TABLE  15 

ACCESSION  USING  QUADRATIC  PROGRAMMING 
Obj. 


Val . 

GURL 

SURF 

SUB 

PILOT 

NFO 

Total 

QP1(Y) 

455.2  | 

I  704 

1404 

554 

1470 

332 

|  4464 

QP3( Y) 

390.6 

|  712 

1416 

563 

1246 

272 

|  4209 

QP  3 ( Y ,X) 

400.8 

|  704 

1380 

556 

1388 

270 

|  4298 

type  and  rank  still  did  well  when  compared  to  the  other  two.  In  fact,  as 
can  be  seen  from  the  actual  accession  values,  all  three  programs  had 
remarkably  close  results.  Since  QP3(Y)  did  the  best  in  terms  of  the  value 
of  the  objective  function,  the  results  from  this  optimization  will  be 
presentee  in  the  remainder  of  this  section. 

While  close,  the  results  from  all  three  programs  have  rather  hign 
penalty  values.  Recall  that  the  minimum  value  that  can  be  achieved 
is  zero,  and  the  actual  magnitude  of  the  penalty  values  shows  that  our 
actual  allocation  is  far  off  the  mark.  Another  way  to  gain  an  apprecia¬ 
tion  of  how  close  the  allocation  is  to  the  targets  is  by  examining  the 
percent  error  compared  to  the  permitted  errors.  Listed  below  in  "able  16 


TABLE  16 


PERCENT  ERROR  3Y  RANK  AND  JOB  TYPE 


1000 

1050 

1110 

1120 

1310 

1320 

1300 

ENS 

351 

0 

19 

54 

125 

-9 

0 

LTJG 

90 

0 

32 

22 

22 

-58 

-12 

LT 

12 

-1 

13 

8 

-18 

-51 

-32 

LCDR 

-26 

-32 

-20 

-46 

-47 

-40 

-44 

CDR 

-14 

-17 

-8 

-35 

-29 

-6 

-19 

CAPT 

14 

31 

24 

23 

0 

0 

47 

are  the  percent  errors  from  the  QP3(Y)  program.  This  table  highlights  a 
problem  that  is  probably  due  to  the  high  rate  of  attrition  of  the  middle 
level  officers.  Notice  that  those  ranks  are  underfilled  (as  indicated  by 
the  minus  sign)  while  the  junior  and  senior  ranks  are  experiencing  an 
overfill.  The  most  serious  problem  exists  in  the  rank  of  Lieutenant 
Commander  and  in  the  job  classes  Submarine  Warfare  Officer  and  Pilot. 

The  model  could  be  made  to  actually  fill  those  shortage  billets  by  moving 
the  permitted  unit  error  in  those  ranks  and  job  classes  closer  to  zero. 
However,  since  the  current  parameters  are  already  closer  to  zero  than  any 
of  the  others,  such  changes  would  only  result  in  worse  overfills  at  the 
high  and  low  end.  A  better  approach  would  be  to  revise  the  sharing  rules 
in  an  attempt  to  get  a  better  allocation. 


VIII.  SUMMARY 


A.  GENERAL 

The  purpose  of  this  thesis  was  to  investigate  and  enhance  the  model 

l 

proposed  by  Grinold  [Ref.  1]  with  special  attention  to  the  optimization 
techniques  and  considerations  involved  with  implementation.  The  model 
looks  at  requirements  and  inventories  of  a  manpower  system  by  classifying 
them  according  to  people  type,  rank,  and  job  type.  By  the  use  of  sharing 
fractions,  which  represent  management  policy,  it  is  possible  to  express 
billets  which  are  usually  defined  by  job  type  as  requirements  for  certain 
types  of  people.  The  same  sharing  rules  also  make  it  possible  to  express 
inventories  of  people  by  the  types  of  jobs  filled  rather  than  by  the 
types  of  people  filling  those  jobs.  Further,  the  assumption  of  steady 
state  makes  it  possible  to  calculate  all  inventories  from  accessions. 

Thus,  billets,  inventories,  and  accessions  are  linked  by  a  unified  and 
rather  elegant  model  structure. 

If  relative  trade-offs  between  different  billets  in  the  form  of 
permitted  shortfalls  or  overfills  can  be  made,  then  it  is  possible  to 
measure  the  error  between  the  desired  billet  structure  and  the  steady 
state  inventories.  Various  optimization  techniques  may  then  be  used  to 
examine  the  long  term  result  of  policy  decisions  regarding  the  allocation 
of  people  to  jobs  and  to  test  alternative  policies. 

This  thesis  has  included  several  areas  of  research  not  presently 
found  in  Grinold  [Ref.  I].  Several  considerations  were  discussed  concern¬ 
ing  the  implementation  of  the  model.  These  include:  level  of  detail  in 
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attribute  definitions,  problems  encountered  in  implementation,  acquisition 
of  model  parameters,  intuitive  explanation  of  weighting  factors,  and 
advantages  of  various  optimization  schemes.  New  programming  formulations 
were  presented  that  used  accessions,  inventories,  sharing  rules,  or 
expected  stage  lengths  as  the  decision  variables.  Elastic  programming 
was  introduced  as  a  generalization  of  the  goal  programming  technique. 
Explicit  solutions  were  derived  for  the  quadratic  programming  model  and 
the  idea  of  additional  equality  constraints  was  proposed.  Other  enhance¬ 
ment  features  included  the  least  squares  approach  to  determine  the 
requirements  by  people  type  and  the  job  sharing  rules  given  the  billet 
structure  and  the  people  sharing  rules,  and  using  least  squares  to 
determine  accessions  from  a  given  inventory. 

In  conclusion,  this  model  allows  the  policy  planner  to  examine  a  wide 
range  of  policy  options  using  only  a  few  simplifying  assumptions  and  a 
modest  amount  of  input  data.  Many  areas  of  policy  planning  may  be 
illuminated  by  simply  placing  a  problem  in  the  framework  of  the  model 
since  this  requires  that  the  analyst  consider  trade-offs  between  various 
people  and  job  types  and  formalize  policy  in  setting  the  values  for  the 
sharing  fractions.  This  ability  to  examine  the  aggregate  effect  of  human 
behavior  on  an  organization's  billet  structure  makes  this  model  a  useful 
tool  for  the  policy  planner. 

B.  AREAS  FOR  FURTHER  STUDY 

There  are  several  areas  related  to  this  model  where  further  research 
should  be  done.  Suggested  areas  include  the  following: 

•  Extend  the  model  structure  to  handle  dynamic  conditions.  If 
the  current  inventory  is  known  and  the  continuation  rates  of 
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personnel  can  be  estimated,  then  a  year-by-year  approach  could 
be  used  to  examine  various  short-range  policy  problems.  This 
would  eliminate  the  limitation  of  steady  state  assumptions  in  the 
current  model . 

•  The  optimizations  rely  heavily  on  reasonable  estimates  of  the 
permitted  percent  error  in  failing  to  fill  all  billets.  More 
research  is  needed  on  how  to  gather,  scale,  and  interpret  data 
that  could  be  used  by  the  model  with  regard  to  these  permitted 
errors. 

•  More  improvement  is  needed  in  adding  cost  considerations  to  the 
model.  The  ability  to  relate  any  analysis  to  "the  bottom  line" 
should  enhance  the  results  of  such  an  analysis.  The  idea  of  a 
multi-attribute  objective  function  that  considers  both  errors  and 
costs  was  presented  in  Section  III.  Another  possibility  would  be 
to  develop  a  means  to  econometrically  express  the  trade-off 
parameters. 

•  The  usefulness  of  the  model  would  be  extended  if  the  results  could 
be  expressed  in  such  a  manner  that  the  impact  of  policy  decisions 
could  be  quickly  evaluated.  One  idea  is  to  use  the  model  to 
create  ratios  or  indices  that  pinpoint  the  criticality  of  various 
allocations.  For  example,  if  the  "optimal"  billet  structure  based 
on  projected  inventories  were  calculated,  then  the  ratio  of  this 
theoretical  structure  and  the  actual  structure  should  provide  a 
measure  of  "how  critical"  certain  billets  will  be  in  the  future. 
Ratios  close  to  one  are  desirable,  and  those  far  away  highlight 
the  oroblem  areas. 
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Since  the  model  is  based  on  steady  state  assumptions,  it  is 
desirable  to  be  able  to  relate  the  results  of  the  model  to  these 
assumptions.  Assuming  certain  functions  that  describe  the  way  in 
which  people  leave  the  organization  could  be  used  to  establish 
such  a  relationship.  The  impact  of  changes  in  continuation  rates 
or  survivor  fractions  on  the  billet  structure  needs  to  be  directly 
related  to  the  sharing  policies  that  will  be  required  to  reduce 
the  allocation  errors,  and  conversely,  changes  in  the  billet 
structure  need  to  be  analyzed  in  light  of  the  behavior  of  people 
in  the  organization. 

Further  model  enhancement  is  needed  in  order  to  introduce  the 
effects  of  stage  substitution  to  the  model.  In  military  organiza¬ 
tions,  such  stage  substitutions  exist  in  the  form  of  "selected  but 
not  yet  promoted."  The  model  has  assumed  that  such  stage  substitu 
tion  (i.e.,  personnel  of  one  rank  filling  a  job  of  a  higher  rank) 
is  negligible  or  has  a  net  effect  of  zero.  The  ability  to  model 
this  would  add  more  credibility  to  the  model  and  subsequent 
analysis. 
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VARIABLE  LIST  FOR  THE  MODEL  AND  PROGRAMS 
Name  Size _ Description _ 


k 


The  index  of  the  manpower  classes.  k=l,2,...,K. 


j 

B 

P 

T 

X 

Z 

A 

F 

G 


The  index  of  ranks  or  career  stages.  i=l,2,...,I. 

The  index  of  the  job  classes.  j=l,2,...,J. 

I,J  The  number  of  people  in  rank  i  required  to  fill  class  j  jobs. 

I,K  The  number  of  required  people  in  rank  i  and  manpower  class  k. 

K, I , J  The  number  of  (i,k)  people  required  in  type  (i,j)  jobs. 

I,J  The  inventory  of  people  in  rank  i  who  are  filling  job  j. 

I,K  The  inventory  of  people  in  rank  i  and  manpower  class  k. 

K,I,J  The  inventory  of  (i,k)  people  who  fill  type  j  jobs. 

K, I , J  The  fraction  of  (i,j)  jobs  that  should  be  performed  by 
manpower  class  k. 

K,I,J  The  fraction  of  type  (i,k)  manpower  who  should  fill  job 
type  j. 


W  I,K  The  expected  length  of  time  a  person  in  manpower  class  k 
wil 1  spend  in  rank  i . 


S  I,K  The  fraction  of  type  k  people  who  will  remain  in  the  organ 

ization  past  rank  i. 

*  I,J  The  permitted  unit  error  for  underfilling  job  class  j 

and  rank  i . 


9  I, J  The  permitted  unit  error  for  overfilling  job  class  j  and 
rank  i . 


5  I ,K  The  permitted  unit  error  for  underfilling  rank  i  and  man 
power  class  k. 


♦  I,K  The  permitted  unit  error  for  overfilling  the  rank  i  and 

manpower  class  k. 

B_  I,J  The  lower  bound  of  permitted  error  for  job  class  j  and 

rank  i . 
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The  upper  pound  of  permitted  error  for  job  class  j  and 
rank  i . 

The  lower  bound  of  permitted  error  for  manpower  class  k  and 
rank  i . 

The  upper  bound  of  permitted  error  for  manpower  class  k  and 
rank  i. 

An  indicator  that  is  1  if  manpower  class  k  can  fill  job  j; 

0  otherwise. 

The  error  in  meeting  (i,j)tn  requirement. 

x.  u 

The  error  in  meeting  (i,k)  n  requirement. 

The  weighting  factor  for  an  (i,j)  error. 

The  weighting  factor  for  an  (i,k)  error. 
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QUADRATIC  FORMS  FOR  THE  MODEL 


1 .  The  General  Form 

The  quadratic  objective  functions  in  this  thesis  have  the  same  qeneral 
quadratic  form  which  can  be  expressed  in  summation  notation  or  matrix 
notation,  if  properly  reindexed.  Consider  the  quadratic  form  with 
var  >les  X(I,1),  A(I,1),  P( 1 , 1 )  and  constant  r: 

1  *iV  - 2 1  »i«i  * r  •  <al> 

i  i 


In  matrix  form  this  becomes: 

X ‘ QX  -  2X 1 P  +  r 

where  Q  is  the  diagonal  matrix: 


Q  = 


r  a  2 

al 


V 

aI  -I 


(a2) 


2 .  The  Special  Form 

If  the  above  general  quadratic  form  has  the  following  characteristics, 
then  a  more  familiar  form  can  be  written  for  it. 


Let  r  =  s  b.  and  p.  =  a^b. 


and  then  (al)  becomes: 


t  aj2xi  -  2=  a(b(x(  +  H>,2  , 


1 


1 


1 


2  (aixi  -  bi'‘ 
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(33) 


In  matrix  notation  this  is: 


X'D'DX  -  2X 1 D ' B  +  B'B 
(DX  -  B)'  (OX  -  B) 
where  D  is  the  diagonal  matrix: 

a. 


’1 


lI  J 


This  special  form  has  several  cases  which  are  of  interest. 
Case  1 

Let  ai  =  1.  Then  (a3)  becomes: 

s  (x.  -  bi )2 
i 

or  in  matrix  form: 

(X  -  B)'  (X  -  B) 

Case  2 


Let  bi  =  a1-ci  .  Then  (a3)  Becomes: 
£  a2  (x.  -  Ci)2 


(a4) 


( a5) 


(a6) 


(a7) 


X'D'DX  -  2X ' D ' DC  -  C'D'DC 
(DX  -  OB) '  (DX  -  DB) 

(X  -  C)'  D'D  (X  -  C)  . 


Case  3 


( a8) 


Let  Dl  =  aiCi  and  x.  =  ifikzk  . 

Since  FZ  =  X,  this  form  is  like  (a7)  and  has  similar  results: 
£  a,2  [(£  fjkzfc)  -  Cij2 

or  in  matrix  notation: 

(FZ  -  C)'  D'D  (FZ  -  C)  . 


(a9) 


,  alO) 
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APPENDIX  C 


APL  FUNCTIONS 


Function  Uses _ Computes 


QUAD  All  variables 


QP3XY  B,  F,  P,  W,  Y,  X,  Z,  A 

e,  <P 


MAKEPG  P,  G  B,  F,  T 

MAKEBF  B,  F  P,  G,  T 

MAKEFP  B,  G  P,  F,  T 

QP3Y  G,  B,  F,  P,  Y,  X,  Z,  A 

W,  8,  * 

INVENTORY  G,  Y,  W  Z,  X,  A 


Syntax  and  Description _ 

QUAD. 

Master  Driver  Program.  This 
function  calls  all  others  through 
the  use  of  menu  options. 

QP3XY. 

Calculates  "optimal"  accessions 
Y,  personnel  inventory  Z,  and 
allocation  A,  given  the  data, 
objective,  and  the  billet  share 
rule  F. 

MAKEPG. 

Reconciles  the  global  variables 
B,  F,  and  T  with  the  global 
variables  P  and  G. 

MAKEBF. 

Reconciles  the  global  variables 
P,  G,  and  T  with  the  global 
variables  3  and  F. 

MAKEFP. 

Reconciles  the  global  variables 
P,  F,  and  T  with  the  global 
variables  B  and  G  using  least 
squares  approach. 

QP3Y. 

Calculates  "optimal"  accessions 
Y,  personnel  inventory  Z,  billet 
staffing  X,  and  allocation  A, 
given  the  data,  and  use  of  the 
people  share  rule. 

INVENTORY. 

Calculates  inventory  Z,  X  and 
allocation  A,  where  accessions, 
people  share  rule,  ana  expected 
stage  lengths  are  given. 


PEOPRCNT  P,  Z 


100(D-Z)/P 


PEOPRCNT. 

Calculates  the  percentage  error 
in  inventory.  Compare  with  *. 


BILPRCNT 

B,  X 

100(8  -  x )/B 

BILPRCNT. 

Calculates  the  percentage  error 
in  meeting  billet  requirements, 
compare  with  e . 

GWAITS 

ALPHA,  S 

W 

GWAITS. 

Given  the  survivor  fractions 

ALPHA  and  stage  definitions  S, 
calculates  W  the  expected  waiting 
time  in  each  stage. 

WAITS 

RET,  S 

w 

WAITS. 

Given  the  retention  rates  RET 
and  stage  definitions  S,  cal' 
culates  W,  the  expected  waiting 
time  in  each  stage. 

ROUND 

ROUND  XXX. 

Takes  any  array  and  rounds 
elements  to  integers. 

LABIJ 

ILAB,  JLAB, 

DP,  CS 

LABIJ  XXX. 

Takes  any  array  of  dimension 
(I,J)  and  labels  rows  and  col¬ 
umns  according  to  ILAB  an  JLAB. 
ILAB  must  have  I  rows  and  JLAB 
must  have  J  rows.  DP  specifies 
number  of  decimal  places  desired. 
CS  specifies  minimum  column 
spacing. 

LABIK 

ILAB,  KLAB, 

DP,  CS 

LABIK  XXX. 

Handles  any  array  of  dimension 
(I,K)  in  the  same  manner  as  LABIJ 

LABKIJ 

ILAB,  KLAB, 
JLAB,  DP,  CS 

LABKIJ  XXX. 

Handles  any  array  of  dimension 
(K,  I,  J)  in  the  same  manner  as 
LABIJ. 

DISPIN 

G,  F,  B,  P, 
ILAB,  KLAB, 
JLAB,  OP,  CS 

DISPIN. 

User  selects  from  a  menu  to 
produce  the  input  variable  with 
appropriate  labels. 

01  SPOUT  X,  Y,  Z,  3, 
ILABS,  KLAB, 
JLAB,  DP,  CS 


DISPOUT 

User  selects  from  a  menu  to 
produce  the  output  variable  with 
appropriate  labels.  Also  displays 
error  arrays  from  3ILPR0NT  and 
PEOPRCNT. 
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NAVY  URL  CONTINUATION  RATES  1 


Years 

GURL 

SURF 

SUB 

PILOT 

NFO 

1 

0.963 

0.971 

0.947 

1.000 

1.000 

2 

0.960 

0.971 

0.938 

1.000 

0.985 

3 

0.930 

0.928 

0.933 

0.990 

0.970 

4 

0.820 

0.809 

0.841 

0.980 

0.919 

5 

0.698 

0.698 

0.697 

0.882 

0.858 

6 

0.816 

0.819 

0.810 

0.692 

0.882 

7 

0.870 

0.860 

0.891 

0.706 

0.859 

8 

0.912 

0.914 

0.908 

0.808 

0.959 

9 

0.947 

0.970 

0.902 

0.818 

0.972 

10 

0.843 

0.833 

0.862 

0.828 

0.934 

11 

0.368 

0.884 

0.835 

0.883 

0.918 

12 

0.922 

0.922 

0.922 

0.933 

0.927 

13 

0.914 

0.926 

0.890 

0.969 

0.971 

14 

0.981 

0.987 

0.970 

0.964 

0.986 

15 

0.974 

0.971 

0.979 

0.954 

0.974 

16 

0.989 

0.990 

0.987 

0.948 

0.974 

17 

0.978 

0.974 

0.987 

0.957 

0.970 

18 

0.956 

0.964 

0.941 

0.950 

0.936 

19 

0.917 

0.901 

0.948 

0.838 

0.954 

20 

0.762 

0.754 

0.778 

0.765 

0.789 

21 

0.331 

0.821 

0.851 

0.774 

0.692 

22 

0.908 

0.883 

0.957 

0.901 

0.914 

23 

0.875 

0.853 

0.920 

0.844 

0.860 

24 

0.911 

0.894 

0.946 

0.852 

1.000 

25 

0.848 

0.847 

0.849 

0.838 

1.000 

26 

0.701 

0.710 

0.682 

0.778 

1.000 

^Source:  OCNO,  OP-Ol. 
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APPENDIX  E 
COMPUTER  PROGRAMS 


VBDISPLAYLU ]7 
7  BDI SPLAY  ;3 


L  1  ] 

* ENTER  SYMBOL  OF  VARIA3LE  TO  BE  DISPLAYED’ 

C  2  ] 

’  TO 

QUIT  TYPE  0* 

l3] 

f  1 

L4] 

’  B 

( BILLET  RE QU I REVERTS)’ 

L  5  ] 

’  F 

( JOB  SHARING  FRACTIOUS)’ 

L  6  ] 

•  T 

( TARGET  ALLOCATIONS)’ 

L  7  j 

’  X 

( CALCULATED  BILLETS  FILLED)’ 

13] 

’  A 

( CALCULATED  ALLOCATIONS)’ 

19] 

’  E 

( PER  CENT  ERROR  IN  MEETING  BILLET  REQUIREMENTS )’ 

.10] 

•3* 

ill]  £0:-(  (S-ft)*’ BFTXAZ’  )/L1,L2,£2.£l,£2,L3 

.12]  -*0 

i.  1  3  ]  L 1 : £ ' LABIJ  ’ ,5 
l14-]  +LQ 

l  1  5  ]  L2i  t’  LABXIJ  •  ,S 

l  1  6  ]  -£0 

t17]  L3  :  z’ LABIJ  ROUND  SILPRCNT 
[1.3]  -*£0 

7 

vsiiPifcwrtaiv 

7  SILPRCNT 

L1.  ]  L  ,431V  PCV//P  1  00 x  (£-3)  f 3  +  3=0 

7 


73J,7Q'P[G]7 
7  Z+BINQP  RIOT 

l  1  ]  3-(  )  .-O.Sxftpl 

L  2  ]  //•*■<¥  ,  [  1  ]  (  p  1  )  ,  0 

l  3  ]  /W((*ffT)  +  .>«tf0*.*,3),[i] 

.  4  ]  -T R 

l  5 ]  INVENTORY 

.6]  Z-*-PF1 

7 


7  3iVS£/h?CVLCI]V 
7  X-BI USE ARCH  P 
Ll]  £PCP:-(FP5>IPCl]-P[2])/0 
i.  2  ]  ,?-(>/P)  +  2 

.3]  -o <3INQP  X)/L) 
l4]  Pl2]-*-X 

,5]  +L00P 

l  6  ]  C  1  :  P  L  1  ]  •*■/ 

L  7  J  -LOOP 
7 
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'/CUANGEZC.1  7 
V  CHaNGE 

l  1  ]  '  *NAT  13  the  name  of  THE  VARIABLE  TO  3E  CHARGED?' 

L  2  ]  A' •►•.!] 

l3  ]  -U*'Q')/SHD 

[4]  -((p,a)=  1  3  5)/C1,C2,C3 

l5J  —ERR2 

16]  C2: tIF  '*/X="PHI_"  '  THE [l  '-L'  ELSE  '-ERR2' 

17]  CS:!/?  '  */X="THETA'  '  '  THEN  '—L  '  '-5^2' 

id]  C*  1  :-<■(  1x  +  /a  =  ’3PFW’  ) /£,Yi?2 

[9]  Li'NHAT  IS  THE  LOCATION  OF  THE  ELEMENT  YOU  NISH  TO  CHANGE?' 
HO]  £0  :■*(  Q  =  x/  (  *  *  p  ,  *  ,X)2XIJ-C)  /ERR4! 
ill]  'INPUT  HEN  '  .X. '  (  *  ,  tlKIJ)  .  '  )  ' 

[12]  INDEX-IKIJZI-I  ] 

L  1  3 ]  £1  :-*(  (p/UV)<r«-H-l) /OUT 

[14]  INDEX-INDEX.'  •,'  .VKIJZI] 

[15]  -£ 1 

[16]  OU?:tX,'L' .INDEX. '1-Q' 

[17]  -END 

[18]  ERR1  INPUT  ERROR.  TRY  AGAIN.' 

[19]  -£0 

[20]  ERR2 : *  SYMBOL  NOT  UNDERSTOOD' 

[21]  'ENTER  AGAIN  OR  TYPE  0  TO  QUIT' 

L  22]  -*1 

[23]  END:' END  OF  CHANGE  PROGRAM' 

7 


HCDIAGlZlN 
7  Z-CDIAG  X'.U'.J 
1 1  ]  IN— pa 
i2]  Z—DG  XZ lJ-l ] 

L  3  ]  £  1  :  *'■(  IV[  2  ]  )  /  O 

l 4 ]  2>Z,[1 ]  *[;«n 

l  5  ]  •*£  1 

7 


7£’«fi  y[d]v 

V  C*P\i 

[1]  'INPUT  TOTAL  NUMBER  OF  RECRUITS  DESIRED' 
L  2 ]  RTOT-Q 

l3]  H-(  (  W1  )-.*UU+.*in)  ,-Q.5*Kpl 

l 4  ]  M,[l](£pl)  .3 

u  5  ]  jf£M((Wl)*.*ff£*.*,3),[1] 

lo]  l-K+YL 

1 
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7  £  6‘  C  Li  ]  7 
7  RR-U'J  X ;  .V 

11]  N*-  pX 

12]  «*i*(  n'  ,N )  p  0 

l  3  ]  Rst'i  1  ;  ]-*-X 

t  **  ]  J  ,  -  i  </-  1 )  $C  1  ]  RR 

7 


VDISPlNiull 

V  DISPIN lS 

Li]  *  ENTER  SYMBOL  OF  VARIABLE  TO  BE  DISPLAYED' 
1 2 ]  '  TO  GO IT  TYPE  O' 

[3]  ’  TO  SEE  ME  HU  TYPE  M' 

l4]  L 4:' 

L  5 ]  *  6  ( BILLET  REQUIREMENTS) ' 

to]  *  P  ( JOB  SHARING  ARRAY)' 

1 7 ]  *  T  ( TARGET  ALLOCATIONS)' 

CU]  ’  G  ( PEOPLE  SHARING  ARRAY)' 

C  9 ]  ’  P  ( PEOPLE  REQUIREMENTS)' 

CIO]  *  7  ( EXPECTED  FURTHER  DURATION)' 

•111]  £0:* 

C 12]  '  • 

L  1 3  ]  *13* 

Cm]  ■*(  1  *p5‘*'i3)  /SRR\ 

it  15]  *-(S='  BFTGPwM'  )  /  L'i  ,L2  ,£5  ,£2  .£6  ,£3  ,£4 

l 1 6 ]  -END 

L 17  3  £  1 :  t ' LABIJ  ROUND  ' ,S 
tl8]  -L 0 

t 19 ]  L2it'LABKIJ  ' ,5 
L  20 ]  -£0 

C  2 1 ]  £3  :  i ' LA3IK  '  ,S 
C  22  ]  -*-£0 

l  2 3 ]  £5 : t  * LABKIJ  ROUND  ' ,5 
C24]  *£0 

t  2  5 ]  £6 : x ' LASIX  ROUND  ' ,S 
l26]  -*£G 

l 27  ]  :  'INCORRECT  ENTRY' 

L28]  END:' END  OF  DISPLAY  PROGRAM' 

V 


VJFCU3V  vr.V£'1*’LlJ]7 

V  WF  4  7  Z9  9-/I9  9  THEN  B 99 

Li]  2*-/i  tlJ  299-,S99Ll+i^99;  ] 

7  7 


7££S£lJ]V 
7  2-4  ELSE  3;F 
1.1]  F*-  (  p  ,  A  )  r  p  ,  3 

l 2 ]  2-(F  ti)  ,  C  J  .  5  ]  Ft/1 

7 
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scv/iirsuil? 

7  u»  AITS  ;A;SF;N;I;AA;33;DD ;  Eli 
'Ll]  H~S 
u 2 ]  3 F^ ALFA 

!_  3  3  N+-1  t  p  3 F 

l4]  I~pa 
l  5  ]  AA+~A  0 . 2  \N 

L  6  ]  3  d-*-.4  «  .  2  *  1  +  \ ,/ 

1 7  ]  33*-d3*~AA 

Id]  DD+^liJ  ,I)pA-lA 
L  9  1  AA"^mA  +  DD*33 

lIO]  .?0<-(7,-7)  +  (7.I+1 )p“l ,1 , 7p0 
L 1 1 3  AA~DD+.*AA 
L  12  3  RR+AA+.xSF 
i  13]  «_—RR 

7 


7i«y£uv:rc'rfrCLl]v 

7  INVESTOR! 

Li]  4W*(7,£)p7 

L2]  7-f/[l]  A-Jx<S)U,l,K)pZ 


VLAdIJLQ]V 

7  Z+-LA3IJ  XuUIJ\I\J\L\HLiMiPP\ND\FS\H\XN\XX'%NS 

[I]  *THIS  PROGRAM  LABELS  All  (7*7)  ARRAY  USING  JLAB  AND  I  LAB  AS 
C  2  ]  *THE  LABELS  FOR  ROWS  AND  COLUMNS  RESPECTIVELY .  DP_  IS  TEE 

L  3 ]  nN UMBER  OF  DECIMAL  PLACES  TO  BE  DISPLAYED  AND  £5  75  THE 
[4]  MINIMUM  SPACING  BETWEEN  COLUMNS. 

L  5  ]  //♦  0 

[6]  I*-{  pX  )C  1  ] 

17]  J—( p X) [ 2 ] 

[a]  ~(I*{oILAB)i\]) /BRR1 

L  9  ]  +U*{pJLA3)l\]) /ERR2 
CIO]  -  L*-{pILAB)L2] 

[II]  HL*-{oJLA3)i2] 

L  1 2  ]  .■!<-(  (  (2,i)p*  *),[!]  ILAB)  .[2](  (7+2)  .2)pt  I' 
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w 27]  2X32:’ ERROR: 
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VIJPUTLl.11 
V  INPUT iX iXIJ itliSiZ 
[1]  'ENTER  THE  NUMBER  OF  MANPOWER  TYPES,  NUMBER  OF  HARKS,  AND ’ 

l  2  ]  '  NUMBER  OF  BILLET  TYPES  AS  A  VECTOR.  ’ 

i  3  ]  XIJ*- u 
[4]  ?  • 
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L 3  J  Li’ 
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[29]  - END 
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L5]  Z-*-\I+-  0 
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